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Abstract 

We consider the situation where a two-level atom is placed in the vicinity of the center of a 
spherical cavity with a large numerical aperture. The vacuum field at the center of the cavity is 
actually equivalent to the one obtained in a microcavity, and both the dissipative and the reactive 
parts of the atom's spontaneous emission are significantly modified. Using an explicit calculation 
of the spatial dependence of the radiative relaxation rate and of the associated level shift, we show 
that for a weakly excitating light field, the atom can be attracted to the center of the cavity by 
vacuum-induced light shifts. 
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I. INTRODUCTION 



Many theoretical and experimental work has been devoted during recent years to the 
so-called "cavity QED" regime, where strong coupling is achieved between a few atoms 
and a field mode contained inside a microwave or optical cavity. In particular, it has been 
demonstrated that the spontaneous emission rate of an atom inside the cavity is different 

0,0,0,0. 



from its value in free space 



This effect can be 



discussed from several different approaches, and here it will be basically attributed to a 
change of the spectral density of the modes of the vacuum electromagnetic field, which is 
due to the cavity resonating structure 0|. This approach is particularly convenient when 
the cavity does not have one single high-finesse mode, but rather many nearly degenerate 
modes, as it is the case in confocal or spherical cavities. More precisely, we will show that a 
"wide aperture" concentric resonator using spherical mirrors with a large numerical aperture, 
can in principle change significantly the spontaneous emission rate of an atom sitting close 
to the cavity center, even with moderate finesse. Similar result were already demonstrated, 
using either a spherical cavity or "hour-glass" modes in a confocal cavity [12|. 

In such experiments, the atom has to sit within the active region volume, which is usually 
of very small size (of order (10A) 3 to (100A) 3 ). In refs 0, 0|, a possible solution was 
implemented by using a narrow atomic beam, and by reducing the cavity finesse in order to 
have an extended area in which a spherical wave is "self-imaged" on itself. However, getting 
large effects will put more severe constraints both on the quality of the cavity and on the 
localisation of the atoms. A different way to implement the proposed scheme in a spherical 
cavity, that we would like to discuss in more detail here, is to use light-induced forces in 
order to attract the atom to the cavity center. A possible implementation could be to couple 
a light beam inside the cavity, and then to use the dipole force to hold the atoms in the 
right position, i.e., close to the cavity center. The effect of strong atom-cavity coupling on 



the dipole force has been studied theoretically [1J, |2JJ , and very interesting effects can be 
expected : since the atomic relaxation will be modified by the cavity, the balance between 
the trapping and heating effects of the dipole trap will be changed with respect to free space, 
which could result in an improvement of the trap itself. 

A good understanding of these effects requires first to know the full space dependence 
of the cavity- induced damping and level shifts. In this letter, we will look at the situation 
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where the atom lies close to the center of a spherical cavity with a large numerical aperture. 
We will show that large changes both in the atom damping rate and in its energy levels can 
be expected, even with a moderate cavity finesse, provided that the atom sits (relatively, 
but not extremely) close to the cavity center [13]. Moreover, we will show that for a weak 



which 



excitating field, the atom can be trapped by vacuum-induced light shifts 
create a force whose spatial dependence is related to the shape of the mode spectral density. 
In the following, we will assume that the cavity damping rate k is much larger than the 
free-space atom damping rate T vac . In that case, the cavity still acts as a continuum with 
respect to the atomic relaxation. This will allow us to treat simply the atom-field coupling 
using frequency- dependant coupling coefficients. For simplicity, we will also assume a weak 
excitation of the atom. This assumption is however not crucial, and possible extensions will 
be discussed at the end of the paper. 



II. LIGHT-INDUCED FORCES IN THE CAVITY QED REGIME 

Let us consider a neutral, slowly moving two-level atom in the dipole approximation. The 
Hamiltonian describing the atom-light coupling is : 

H = H f + H a - D.E(r) (1) 

where D is the atom dipole moment and E is the electric field. The free field and atom 
hamiltonians are respectively Hf = Y,k + 1/2) and H a = p 2 / (2m) + %uj S^ . Here, r 

and p are the atom's position and momentum, m is its mass, and the dipole operator can be 
written D = d(e lulLt + e~ lu}Lt S^) , where and are the usual two-level rising and 
lowering operators in the frame rotating at the angular frequency ojl of an externally applied 
laser (we have = (S^S^ — S^) /2). The electric field operator is expanded as 
usual on a basis of orthogonal modes, but we do not specify that these modes should be 
plane waves. One has therefore : 

E(r, t) = i 52 e k(r) a k (t) e*** % + h.c. (2) 
k 

where k is a mode label which include polarization, and e&(r) is the contribution of mode 
k to the field at point r close to the cavity center. It will be convenient to refer optical 
frequencies to a reference ujl, which can be the frequency of an externally applied laser as 
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said above, and to define : 

A(r,t)=i£e fc (r)a fc (t)e-^-^ (3) 

k 

so that : 

E(r, t) = A(r, t)e~ iuJLt + A f (r, t)e iuJLt . (4) 

Moving to a frame rotating at the frequency ujl and using the rotating wave approxima- 
tion, the interaction part of the Hamiltonian can be written : 

H mt = -S (+) d.A - d.A^ (5) 

From the expression of the hamiltonian one can simply get the force acting on the atomQ] : 

F(r, t) = ^ = -VH = V{S i+ \t) d.A(r, t)) + h.c. (6) 

In this expression, the only space-dependant part is efe(r), so the force can also be written : 

F(r,t) = i S^ + \t) V[d.e fe (r)] a k (t) e^™)* + h .c. (7) 

k 

From the Hamiltonian, one gets the Heisenberg equations for the field operators ai. jl7|. 
which can be integrated formally, yielding (d.e k is taken real) : 

a k (t) = a k (t ) + ^ I''' drS^it - r ) e <to-^)(*-r) (8) 



h Jo 

Using this result in the definition of the fields operators, one sees that the field splits in two 
parts E = E + E s , where E Q and E s are the well-known "vacuum field" and "source field" 
terms, which do not commute \^\ . Similarly, the expression of the force splits in two parts : 

F = F G + F s (9) 

By choosing the normal ordering when separating the two non-commuting terms, the first 

part F gives the usual expression of the light-induced force [l^. The second term F s is 

zero in the absence of a cavity, because the gradient of the source field is zero at the dipole 

place. However, as we will show now, this is no longer true inside a cavity. Still assuming 

normal ordering, one gets : 

zV[d.e fc (r)] d.e fc (r) p ^ (+) _ + ^ 

u n Jo 



F„ ^ 

k 

E i (d.e k (r)f p drS{+){t)s ^ {t _ r) e _ i( . fe _. i)7 

I Zilb J 



V 



+ h.c. (10) 



In order to get the physical meaning of this integral, let us write the Heisenberg equation 
for and use the expression of in order to get : 

^ = i5 L S^ - 2tfW^ + 2 £ f~ t0 drS^(t)S^(t - r)e<^>, (11) 

where ^ = lol — uo . The integral appearing here has the same structure as the one appearing 
in eq. blunder the gradient, and corresponds to the well known relaxation and light shift 
terms in the evolution of the dipole components. It can be calculated as usual using a 
Markov approximation |l7j], which is possible here because as said above we assumed that the 
cavity features are wide compared to the free-space relaxation rate of the atom. In the weak 
excitation limit, which is assumed here, one can approximate the evolution of the correlation 
functions by their free evolution during the short memory time of the reservoir. Note that 
an extension should be made to the strong excitation regime by considering relaxation in a 
dressed state basis 201 . Making the approximation S^~\t — r) = S^~'(t) e^ UL ^^ T during 
the correlation time [17 1, one can define the quantities : 

k h Jo k h V \U k -LU JJ 

r(r)/2 + zA(r) (12) 



where one has therefore : 

9-7T 

r(') = jEN')f«i-».) (is) 

A W =£(^ P (_L-) (14, 
Using these definitions and S^S^ = -S^/2, eq. EU becomes as expected 

^ = (i5 L - iA(r) - ^) fiH - 2tf«^ (15) 

so that the atomic frequency ujq becomes c^o + A(r). We note that the free-space value of 
A(r) is a diverging quantity, which is usually assumed to be absorbed in the definition of 
the atomic levels; therefore, one considers here only the (finite) change of A(r) with respect 
to this free-space value, that will be denoted A'(r) : 

A'(r) = A cat) (r)-A mc (r) (16) 
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In the expression of the extra term in the force, one can use either A'(r) or A cav (r), since 
the gradient of A mc (r) is zero anyway, and one obtains finally : 

F s = h/2 (S^S^ V(-A / (r)+?r(r)/2) + V(-A / (r)-2r(r)/2)S (+) 1 S ( - ) ) 

= -V(ftA'(r)) n e (17) 

where IT e = S^S^ = + 1/2 is the excited state population. Therefore the gradient 
of the cavity-induced light shift creates a force, which can attract the atom towards the 
cavity center for an appropriate choice of the atom-cavity detuning. Physically, it makes 
sense that F s is also proportionnal to the excited level population. In order to characterize 
more precisely the behaviour of this force, one needs now the explicit space dependence of 
A'(r). The same calculation will also give the space- dependent damping T(r)/2, which will 
be useful for calculating the steady state and evolution of the atomic operators. 



III. CAVITY-INDUCED RELAXATION RATES AND LIGHT SHIFTS 

For definiteness, we will consider the case of a spherical cavity of radius R and of reflec- 
tivity and transmittivity coefficients p and r, with p 2 + r 2 = 1, and r 2 = T. We will assume 
that kR ^> 1 (typically kR = 10 5 ), and a moderate cavity finesse (in the range 10-100). A 
crucial parameter is the solid angle subtended by the cavity, that will be denoted AQ cav . 
For instance, a cavity half-aperture angle of 45 degrees gives A ^ Q " = 0.3, and therefore 
A ^ ac = 0.7 as the fraction of space still occupied by vacuum modes. All parameters quoted 
above seem accessible from an experimental point of view, and we will show now that they 
allow one to get quite significant cavity-induced effects. 

In order to calculate the explicit space and frequency dependence of the relaxation rate 
and level shift, we have followed two parallel approaches, which are described in detail in 
another publication [13[ . The first one is to solve explicitly the field equations in a spherical 
geometry, taking into account the considerable simplifications which appear since we are 
interested in the field at distances smaller than ~ 100A off the origin. In geometric optics 
this involves light-rays with an impact parameter smaller than 100A, and therefore, in the 
multipole expansion of the field, it will be sufficient to consider harmonics up to order 
/ < 100. Actually, up to 300 harmonics have be used in order to check consistency. Another 
very important point is that the continuity equations of the fields on the mirrors must be 
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expressed at a distance R ~ 1 cm. With A = 780 nm, this corresponds to kR ~ 80000, which 
is quite large, and allows one to use asymptotic forms of the solutions on the mirrors. Under 
these assumptions, it can be shown that the expression of A'(r) and T(r)/2 can be given an 
explicit operatorial form in the space of mode functions, and then evaluated numerically in 
a spherical harmonics basis Q|. 

Besides this "exact" calculation, we have also looked for an approximate solution, inspired 
by ray-optics considerations, and eventually checked by comparison with the complete nu- 
merical calculation. Using this cross-checking method, we obtain finally that the effect of 
the cavity can be described to a very good approximation by the following formulas: 

r(r) f dQ 3_ ,d.Q. 2 ., T 2/lA T . 2 * 

— ^ = / — ~(l-(— r-) )(< ^ttt7 cos 2 (kQ.r) + — — sin 2 (M2.r)) (18) 

A'(r) f dQ 3. ,d.f2. 9 . . p sin(20) 9/ ,a n /? sin(20) 9/ ,a „ / s 

—A^= / (l-(__) 2 )(Jl \Z> C os 2 (kQ.r) - sin 2 (kQ.r)) 19 

T mc J A* 4tt 2 v v d ' pe 2 ^\ 2 v ; |l + pe 2 ^| 2 v >) \ ) 

where the notation Q describes a direction in space, while is a cavity detuning parameter 

that will be detailed below. These expressions have a straighforward interpretation, because 

they appear basically as integrals over the direction of light rays : in the integral over 

the directions, p is the mirror reflectivity for rays subtended by the cavity, and is zero for 

rays outside the cavity solid angle. The different factors appearing in the integrals can be 

interpreted in the following way : 

• The first factor under the integral corresponds to polarisation effects, taking into 
account the transverse character of the field. 

• The second (resonance) factor is of the usual Fabry-Perot form, where <fi is the cavity 
phase shift which includes first a term O = ujqR/c. The complete calculation shows 
that, in order to obtain a correct result outside the cavity center, <ft should include 
also a contribution from spherical aberrations, that is : = 0o + ~ ^ • This 
second term corresponds to the extra phase shift experienced by rays going through 
point r while propagating along the Q direction. The resonance factor has different 
expressions for the damping and the lamb shift terms, which correspond respectively 
to the active and reactive parts of the coupling. This is clearly apparent from the 
integrals of eq. EI] and El which involve either a delta function or a principal part. 
In the first case, the integration is trivial, and yields the resonance term of eq. [HE 
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while in the second case the result is obtained by contour integration jl^, and gives 
the (dispersive) second term of eq. [Till 

• The third term under the integrals is the stationnary wave pattern corresponding either 
to odd modes (with an anti-node in the center and a cos 2 (Ml. r) space dependence) or 
to even modes (with a node in the center and a sin 2 (/cfi.r) space dependence). 

• Finally, the integration over the mirrors is conveniently performed in spherical coor- 
dinates, by taking the z axis along the cavity axis, and varying the azimuthal angle 
9 from to 9 mirror = 9 m . Improved accuracy (better than 1%) is obtained if one 
takes into account the fact that the rays which would be reflected near the edge of the 
mirror are actually lost due to diffraction and fail to do as many round-trips as the 
other ones. We have shown [13] that this effect can be taken into account very simply 
by decreasing 9 m to 9 e ff = 9 m — 59, with 59 = 1 / V kRT for symmetrical mirrors. 

The first results which can be obtained from the previous formulas are obviously the shift 
and damping at the cavity center, as a function of the atom-cavity detuning. For a dipole 
orientation parallel to the cavity axis, we obtain straightforwardly : 

T par (0) _ Att vac n | sm 2 9 m ^ | AOg^ cos 9 m ( 1 ~l~ cos 9 m ) T 

A par(°) _ ^cav cos 9 m (1 -j- cos 9 m ) p sin(20 o ) 
T vac 4tt 1 2 j |i-pe 2 ^|2 I ) 

while for a dipole orientation perpendicular to the cavity axis, we have : 

T perp (0) _ AQ vac sin 2 ^ AQ cav cos6> m (l + cos6> m ) T 
T vac 4vr 1 4 ' 4vr 1 4 ' \ 1 - pe 2i ^ |2 1 > 

A perp(Q) _ Attcav n , cos 9 m ( 1 + cos 9 m ) p sin(20 o ) 

T vac - An {l+ 4 } \l-pe^\2 ^ 

We note that these expressions yield for a randomly oriented dipole : 

T av (0) An vac A£l cav T A' av (0) An cav p sin(20 o ) 



(24) 



T vac ' 4vr 4?r |l-pe 2 ^o| 2 ' T vac 4tt \l-pe 2i ^\ 2 

which have a straightforward interpretation in terms of resonant enhancement of the rays 
subtended by the cavity. We note that these results are the same as those given in ref. up 
to factor two resulting from the fact that this reference was considering spatially averaged 
values rather than the peak value at the cavity center (see below for the space dependence). 
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These functions are plotted on fig. Q]fbr ^ = 0.3 and p = 0.98. It can be seen that very 
significant effects occur for these quite reasonable parameters, yielding more than 30-fold 
increase in the damping rate at the cavity center. 



FIG. 1: Normalized damping r(0)/T mc (left) and level shift A'(0)/r mc (right) at the cavity 
center, as a function of the atom-cavity detuning (u>q — co cav ) normalized to the cavity linewidth. 
The amplitude reflection coefficient of the mirrors is taken to be p = 0.98, and the numerical 
aperture of the cavity is 0.7. The upper curves correspond to a dipole oriented perpendicular to 
the cavity axis, and the lower curves to a dipole oriented along the cavity axis. 

We can then look at the results as a function of space for a given frequency, which are 
essential for the present paper. Two atom-cavity detunings are specially worth looking at : 
the resonant frequency at the cavity center, which yields maximum change in the damping 
rate but no cavity shift, and frequencies detuned by plus or minus half a cavity linewidth, 
which yield maximum cavity shifts. As an example, the results for the damping rates are 
plotted on fig. |21 The general values of the damping and level shifts for arbitrary values of 
the reflectivities px, p2 of the two mirrors are given in jl3j |. Appendix C. 

IV. DISCUSSION 

From the results of the previous section one can deduce some features of the "vacuum- 
induced" force given by eq. El For simplicity, we consider an atom at point r with zero 
velocity. The simplest configuration in which this force should be dominant is a weak 
stationnary wave resonant on the atom, but detuned from the cavity. In that case, both the 
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FIG. 2: Spatial variation of the atomic damping r(r)/r mc on the cavity axis, obtained from the 
analytical formula given in the text. The atom-cavity detuning is taken equal to zero at the cavity 
center (see Fig. and the horizontal axis unit is 1/k. The amplitude reflection coefficient of 
the mirrors is p = 0.98, and the numerical aperture of the cavity is 0.7. The upper left curve 
corresponds to a dipole oriented perpendicular to the cavity axis, the upper right curve to a dipole 
oriented along the cavity axis, and the lower curve is a zoom close to the cavity center. 

usual scattering and dipole force are zero, while the cavity-induced light shift is maximum 
for a detuning of about half a cavity width. It can be seen easily that the force is attractive 
when the cavity is detuned to the blue side of the atomic resonance (u>o < oj cav ) . This can 
be understood physically by realizing that the cavity "repels" the excited state level, and 
that this level has to go down in order to get an attractive potential. The corresponding 
potential wells are a few hY (see eq. ITTj) . which is quite significant if cold atoms are used. 

Besides the effect described here in the weak excitation regime, another possibility is to 
couple a laser beam inside the cavity, with a red atom-laser detuning, in order to create 
a dipole trap attracting the atom towards the cavity center. A very interesting effect can 
then be expected : since the atomic relaxation will be modified by the cavity, the balance 
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between the trapping and heating effects of the dipole trap will be changed with respect to 
free space, enabling new regimes where the dipole force can be both attractive and damped. 
In that case, the calculation above should be extended to moving atoms (non-zero velocity) 
in order to calculate friction forces, as well as to calculate the diffusion coefficient. For very 
cold atoms, the quantum aspect of motion should also be included; this could be done in 
principle since we have obtained the explicit space dependence of the trapping potential. 



V. CONCLUSION 



As a conclusion, we have shown that macroscopic cavities with large numerical apertures 
are interesting candidates for cavity QED experiments in the optical domain. The possibility 
to use light-induced force to hold the atom close to the cavity center is quite attractive, and 
it offers the possibility to study a well defined quantum system, including its external degrees 
of freedom. The expressions that we have obtained are general, and by using the results 
of ref. 3] they can be applied to any kind of low-finesse cavit y, i ncluding the case of a 



"half-cavity" with only one mirror imaging the dipole onto itself 
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